EPJ manuscript No. 

(will be inserted by the editor) 



Non-Markovian entanglement dynamics in coupled 
superconducting qubit systems 

. . . Wei Cui 1 * 2 , Zai-Rong Xi 1 a , and Yu Pan 1 ' 2 
O ' 

f*"^ , 1 Key Laboratory of Systems and Control, Institute of Systems Science, Academy of Mathematics and Systems Science, Chinese 

(N 

5^ ' Academy of Sciences, Beijing 100190, P. R. China 

2 Graduate University of Chinese Academy of Sciences, Beijing 100039, P. R. China 

Received: date / Revised version: date 

Abstract. We theoretically analyze the entanglement generation and dynamics by coupled Josephson junc- 

-I— > ■ 

tion qubits. Considering a current-biased Josephson junction (CBJJ), we generate maximally entangled 



3 . 



states. In particular, the entanglement dynamics is considered as a function of the decoherence parameters, 
such as the temperature, the ratio r = uj c /u)q between the reservoir cutoff frequency ui c and the system 



oscillator frequency 6Jq, and the energy levels split of the superconducting circuits in the non-Markovian 

m ; 

master equation. We analyzed the entanglement sudden death (ESD) and entanglement sudden birth 

m 

(ESB) by the non-Markovian master equation. Furthermore, we find that the larger the ratio r and the 
thermal energy fcsT, the shorter the decoherence. In this superconducting qubit system we find that the 

o' 

entanglement can be controlled and the ESD time can be prolonged by adjusting the temperature and the 
superconducting phases $ & which split the energy levels. 



PACS. 03.67.Mn Entanglement measures, witnesses, and other characterizations - 85.25.-j Supercon- 
ducting devices - 42.50.Lc Quantum fluctuations, quantum noise, and quantum jumps 

1 Introduction classically, the properties of one particle can depend on 

those of another (typically distant) particle in a way that 

Entanglement is one of the remarkable features of quan- 

only quantum mechanics can explain. Einstein-Podolsky- 

tum mechanics. Briefly, entanglement refers to correlated 

Rosen (EPR) entangled states, probably the simplest and 

behavior of two or more particles that cannot be described 

a n " ~7, r„. most interesting entangled states, have been employed not 

Present address: zrxi@iss.ac.cn ° ° J 



2 Wei Cui et al.: Non-Markovian entanglement dynamics in coupled superconducting qubit systems 

only to test Bell's inequality, but also to realize quan- two-qubit coupling and a controllable coupling mechanism 
turn cryptography, quantum teleportation, and quantum have been realized [19120) . However, due to the difficulty 
computation [1 2 3 1415] . It has become clear that entangle- to decouple the qubits from the environments in solid- 
ment is a new resource for tasks that cannot be performed state systems, only the two-qubit entanglement of the su- 
by means of classical resources |6I7) . Although many re- perconducting qubits has been observed in experiments 
suits have been obtained (for example, see the review pa- [3T] , while the entangled states up to eight [35] or six 
pers [819] ). the theory of quantum entanglement has open photonic [23] qubits have been experimentally reported, 
questions, like (i) how to optimally detect entanglement Therefore, the future generation of multi-particle entan- 
theoretically and practically; (ii) how to reverse the in- gled states for solid state quantum computation would 
evitable process of degradation of entanglement; and (iii) be a significant, step towards quantum information pro- 
how to characterize, quantify and control entanglement, cessing. A challenge is quantum decoherence, because any 
The focus of this paper is a theory study of entanglement pure quantum state used evolves into a mixed state due 
dynamics and control in a coupled Josephson junction sys- to the unavoidable interactions with the environment. De- 
tem. coherence describes the environment-induced suppression 



Superconducting quantum circuits jlO] are the sub- 
ject of intense research at present. Josephson devices can 
serve as quantum bits (qubits) in quantum information 
and that quantum logic operations could be performed 
by controlling gate voltages or magnetic fields 10 . More- 
over, for its scalable and macroscopic property, Joseph- 
son junctions offer one of the most promising candidate 
served as hardware implementation of quantum computers 
[11I12I13I14I15I16I17) . The charge, flux, and phase qubits 
are three basic types of superconducting qubits [TU] de- 
pending on which dynamical variable is most well de- 
fined. Operations with multiple superconducting qubits 
have also been performed. Several types of Josephson- 
junction qubits have been proposed and explored in the 
laboratory. The first solid-state quantum gate has been 
demonstrated with charge qubits [18[ . For flux qubits, 



of the quantum mechanical coherence properties and in- 
terference ability, which transforms the quantum system 
into classical one. The description of this process requires 
us to take into account not only the degrees of freedom 
of the system of interest, but also those of the environ- 
ment. Decoherence of the Josephson junction qubits is 
considered to be the major impediment for quantum logic 
gate operations. Thus, short coherence times limit both 
the manipulation of the qubit state and information stor- 
age. In all superconducting qubits, both the spectrum of 
charge noise and the critical current fluctuations as it is 
display a 1/ f behavior at low frequencies. Moreover, both 
charge noise and critical current fluctuations can be phe- 
nomenologically explained by modeling the environment 
as a collection of discrete bistable fluctuators, representing 
charged impurities hopping between different locations in 
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the substrate or in the tunnel barrier. Many efforts [24] tions in the Lindblad form are characterized by the fact 

searched for decoherence mechanisms that suppress the that the dynamical group of the system satisfies both the 

decoherence time. semigroup property and the complete positivity condi- 

/-v -v. -v. , j tion. However, in the superconducting qubit system the 

Un trie path to quantum computing, superconduct- ° J 

i -i rrrn i i ,, , . . Markovian approximation is not justified 131 321. In re- 

mg qubits [lUj are clearly among the most promising can- J ' — — 1 

j., , n c i ap i a^ i a^ i i j xt„ x u i i i cent years, non-Markovian quantum dissipative systems 

didatcs. References 25 26 27 showed that such two-level j j m r j 

■ V, i i, a u-x ii t I33I34I35I36I37I38I39I40I41T421 have attracted much atten- 

systems can themselves be used as qubits, allowing for 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

ln n j . . i r tion due to its fundamental importance in quantum infor- 

a well controlled initialization, universal sets of quantum 1 ^ 

gates, and readout. Thus, a single current-biased Joseph- ma ^ 10n processing . 

son junction can be considered as a multi-qubit register. In this paper, we will focus on the dynamics of the gen- 
It can be coupled to other junctions to allow the appli- erated superconducting entangled state in a non-Markovian 
cation of quantum gates to an arbitrary pair of qubits environment by using the master equation method. It dc- 
in the system. These results [251261271 indicate an alter- termines how much quantum information can be reliably 
native way to control qubits coupled to naturally formed transmitted over the noisy quantum channels. We then 
quantum two-level systems, for improved superconducting discuss schemes for controlling the entanglement between 
quantum information processing. Nevertheless, the path is the superconducting quantum circuits. Nowadays, quantum- 
long, and there are quantitative technological obstacles to state engineering, i.e., active control over the coherent dy- 
be overcome, notably increasing the decoherence time, im- namics of suitable quantum-mechanical systems, has be- 
proving the fidelity of the read-out, and strengthening the come a fascinating prospect of modern physics [25 26127) . 
entanglement distillation in its dynamics. Quantum decoherence and entanglement control pave the 

,„•,,. xu xu r x x | AA | AA | AA i way for future long coherent time quantum information 

Within the theory of open quantum system 28 29 30 J ° H 

j. . . j i j -u j t, xi. processing and computation, 

the dissipative dynamics can be described by the mas- 

ter equation of the reduced density matrix. In general This paper is organized as follows. We first introduce 

this equation is obtained by tracing over the environ- superconducting qubits decoherence and the quantum non- 

ment variables, after performing a series of approxima- Markovian master equation for driven open quantum sys- 

tions. The Born-Markovian approximation is usually used terns. In Sec. II, we shall briefly recall the physics of the 

in deducing the master equation, which neglects the corre- Cooper pair box, and give the model of the superconduct- 

lations between the system and the reservoir. The Marko- ing qubits interacting with the bath. In Sec. Ill, we con- 

vian approximation leads to a master equation which can sider the entanglement dynamics by the non-Markovian 

be cast into the so called Lindblad form. Master equa- master equation. Both entanglement sudden death (ESD) 
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Fig. 1. Color online) Schematic diagram of the coupled-qubit circuit with a biased-current source of impedance Ziyi). Two 
Josephson charge qubits are controllably coupled to a common current-biased Josephson junction, which operates as a Josephson 
phase qubit and acts as a coupler. Two Josephson qubits and current-biased Josephson junction with the electromagnetic 
environment represented by the impedance Z(iS). 



and entanglement sudden birth (ESB) are analyzed by 
numerical simulation. Conclusions are given in Sec. IV. 

2 The model 

In this work, we present an experimentally implementable 
method to couple two Josephson charge qubits and to gen- 
erate, detect and control of macroscopic quantum entan- 
gled states in this charge-qubit system. Let us study the 
superconducting circuit shown in Fig.l, where two sin- 
gle Cooper pair boxes (CPBs) are connected via a com- 
mon bus, i.e., a current-biased Josephson junction (CBJJ). 
Each qubit consists of a gate electrode of capacitance C g 
and a single Cooper pair box with two ultrasmall identi- 
cal Josephson junctions of capacitance Cj and Josephson 
energy Ej, forming a superconducting quantum interfer- 
ence device (SQUID) ring threaded by a flux <P and with 
a gate voltage V. The superconducting phase difference 
across the fcth qubit is represented by k = 1,2. The 



large CBJJ has capacitance Cf,, phase drop 8b, Joseph- 
son energy Eb, and a bias current lb- The reason for this 
choice is that this circuit can be easily generalized to in- 
clude more qubits, coupled by a common CBJJ. For the 
detail of controllable coupling of superconducting qubits, 
one can refer to |43l44l45l46l47l48l49l50l5I] .The Hamilto- 
nian without dissipation can be written as 

2 

H s =Y J ^E c {h k ~n gk ) 2 -Ej{$ k )cos6 k + H kb + Hb, (1) 
fc=i 

where 

H kb = TrC g Ej($ k )0 b sin4/ [$o{C g + 2C 7 )] , 
H b = Ql/(2C b ) - E Jb cos 6 b - $ I b /{2n9 b ). 

Here E c — e 2 /[2(C g + 2Cj)] is the single-electron charging 
energy of a single CPB. For simplicity, we assume that E c 
and Ej are the same for the two CPBs. C g and C g are the 
capacitances of the gate electrode and identical Josephson 
junctions, respectively. The Josephson energy Ej(<P k ) of 
the fcth dc SQUID is Ej(<P k ) = 2E°jCOs{<K$ k /$ ), where 
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Ej represents the Josephson energy of a single Joseph- standard Pauli operator and uib = Eb\ — Ebo the eigenfre- 
son junction, <P k denotes the external flux piercing the quency. 

SQUID loop of the j-th CPB, and <Pq is the flux quantum. Under the rotating-wave approximation, the Hamilto- 
Qb = 2-Kpb/^o is the operator of charges on the CBJJ. Ejb nian H can be rewritten as 

and C b = Cjb + ZUiC^+C;, 1 )- 1 = C Jb + 2{{2C^ + 2 f E Jk (<P k ) {k) _ E ck (n gk ) (fc A (b) 

■"s — 2_^i 1 2 U z 2 U x I UJb<Jz 

Cg 1 )^ 1 are the Josephson energy and effective capaci- k=l 

+J {a\a 2 _ + a%a x _) , (2) 

tance of the CBJJ, respectively. Also, <Pq = h/(2e) is the 



flux quantum. The operators h k and § k satisfy the com- with 
mutation relations [0 k ,h k ] = i (we assume % = 1), which 
describe the excess number of Cooper pairs and the effec- 
tive phase across the junctions in the fc-th CPB, respec- 
tively. In addition, the phase operator 8b for the CBJJ 



ai k) = |+M-| + 

4 k) = -i\+)kk(-\+i\-)kk{ 
ai k) = \+) kk (+\-\-) kk (-\, 



(3) 



anda£ fc) = (a™ ± ia^ /2, where |+) fc = (|0> fc +|l) fc )/V2 



. and \-) k = (|0) fc - |l))/-\/2. Note that, the charging en- 

and its conjugate pb satisfy another commutation relation 

- . ergy E ck (n gk ) = 2e 2 (f - 2n gk )/C k of the fc-th Joseph- 

[Vb,Pb\ = i- 

son charge qubit can be switched off by setting the gate 

Suppose that the CPBs are biased at the charge de- TA , . , , n , j- 

° voltage V k such that n gk — 1/2. Also, by adjusting the 

generate point, such that n gk = C gk V k /(2e) = 1/2 (when flux ^ ^ Josephson energy of the fc _ th qubit 

V k = e/C gk ). The two energy levels of the fc-th CPB = ^ cos( ^ fcM)) can be set to the strongest 

corresponding to n k = 0,1 are close to each other and „ N , ,, , ,. ^ . 

^ 6 ' coupling (<£ fc = 0) and the decoupling (<P k = $0/2). This 

far separated from other high-energy levels. In this case, n . in , r ,, . 

00 achieves the controllability of the present quantum circuit. 

they behave as effective two level systems (with the basis „> . , T . ,, , , ,, 

v I he cocfhcient J is the coupling strength between the fc-th 

i\ tk) = \n k = 0), I Ik) = \n k = I)|) [H]. It is well known , ., , T . ... T . . ,„ . 

u 1 ' 1 1 nil — 1 qubit and Josephson junction with J = \i\ 2 /E c , \ k = 

that the CBJJ can be approximated as a harmonic oscil- „ ^ n • i^. m\ 
^ 2Ej sin(<P k /2). 

lator 152 1151531541 . if it is biased as h <C In — 2nEjh/ ( l > n. . , . r , , 

1 — — — — 1 Obviously, it we assume the two qubits are both at 

Here, we consider a very different case, i.e., the biased dc ,, , n . , , , , _ /ri 

the charge degenerate point n gk = 1/2, and = <Pn/2 

current is slightly smaller than the critical current 7n, . , . ... , , N , N , > ,-, , ,. 

" ° J u ' with the separated initial state |<?)fc|g)i|e)2, the evolution 

and thus the CBJJ has only a few bound states. The two r ,, n . , .. . . 

operator of tne corresponding two qubit system is given 

lowest energy states |0&) and \lb) are selected to define a ^ 
Josephson phase qubit acting as a two-level date bus. Un- 

U(t) = exp \~iJt{(T 1 + <7 2 _ + a\a x _) - iu} b ta^] , (4) 
der such condition, the Hamiltonian of the CBJJ reduces L -I 

to Hb = TioJbcrl: with = | Ob) (Of, | — |lf,)(lf,| being the and the system's evolution is given by 
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U{t)\9)b\9h\e)2 = e-^? ) \g) i e- iJt (<' 1 +°'- + ° a +° X -'>\g) 1 \e) 



turn dynamical semigroup. The Hamiltonian Hs(t) de- 
scribes the coherent part of the time evolution. Non- negative 

2 

quantities 7 m play the role of relaxation rates for the dif- 

= (cosujbt + i sin LUi,t)\g) i, [cos Jt\g)i\e)2 — i sin Jt\e)\ |flf)a] ■ 

ferent decay modes of the open system. The operators C m 

( 5 ) 

are usually referred to as Lindblad operators which rep- 
If we choose Jt = tt/4, we can obtain the maximally en- resent the various decay modes, and the corresponding 
tangled state \&) = {\g)i\e) 2 - i\e)\ \g) 2 )/V2. density matrix equation © is called the Lindblad mas- 

ter equation. The solution of Eq. © can be written in 
3 Entanglement dynamics terms of a linear map V(t) = exp(Ct) that transforms the 

initial state p(0) into the state pit) = F(£)p(0) at time 

3.1 Markovian and non-Markovian master equations 

t. The physical interpretation of this map V(t) requires 
We account for the dissipation due to electromagnetic flue- that it preserves the trace and the pos itivity of the den- 
tuations. They can be modeled by an effective impedance gity matrix ^ The most important physical assumption 
Z(w), placed in series with the voltage source and pro- which underlies Eq Q is the validity of the Markovian 
during a fluctuating voltage. The impedance is embedded ap p roximation f short environmental correlation times, 
in the circuit shown in Fig.l, which further modifies the with tMs a p proximation) the environment acts as a sink 
spectrum of voltage fluctuations. In general, the environ- for the system informa t ion . Due to the system-reservoir 
ment consists of a large set of harmonic oscillators, each inter action, the system of interest loses information on its 
of which interacts weakly with the system of interest, i.e. gtate into the environment, and this lost information does 
Hbath = J2i ^ T^ x i ■ ^ S mentioned in Sec. 1, we not p i ay any f ur ther role in the system dynamics. 



will use the master equation method to study the entan- 
glement dynamics and control. 

The analysis of the time evolution of open quantum 
system plays an important role in many applications of 
modern physics. With the Born-Markovian approxima- 
tion the dynamics is governed by a master equation of 
relatively simple form |30j 



If the environment has a non-trivial structure, then the 
seemingly lost information can return to the system at a 
later time leading to non-Markovian dynamics with mem- 
ory. This memory effect is the essence of non-Markovian 
dynamics [33I35I36I37I38I39I40I41142] . which is character- 
ized by pronounced memory effects, finite revival times 
and non-exponential relaxation and decoherence. Non-Markovian 



^p{t) = -i[H s itlp{t)]+Y,lrnV[C m ]p{t), (6) . . 

at — ' dynamics plays an important role m many fields ot pnysics, 

with a time-independent generator in the Lindblad form, such as quantum optics, quantum information, quantum 

This is the most general form for the generator of a quan- chemistry process, especially in solid state physics |41) . As 
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a consequence the theoretical treatment of non-Markovian son junction, the bipartite dynamics is 



quantum dynamics is extremely demanding. However, in 



dp(t) 



H s ,p}+J2(r 1 D[a^}p + r 2 D[a+}p), (8) 



order to take into account quantum memory effects, an ^ k=i 

integro-differential equation is needed which has complex where A = A(t) + 7(f) and T 2 = A(t) - j(t). The time 

mathematical structure, thus preventing generally to solve dependent coefficients A(t) and are diffusive term and 

the dynamics of the system of interest. An appropriate damping term, which can be written as follows 
scheme is the time-covolutionless (TCL) projection oper- 
ator technique [30141142] which leads to a time-local first 

order differential equation for the density matrix. j(t) = / dr/i(r) sin (wt), (10) 

Jo 

By tracing out the bath degrees of freedom, we find 
for p a non-Markovian evolution equation 



t 

A(t) = I drk(T)cos(ujT), (9) 







with 



—p{t) = -i[H s (t), p(t)\ + A m {t)V[C m (t)]p(t). (7) 



with the super-operator V\L\p = LpL^ - \L^Lp - \pl^L. ^(r) =2 dwJ(u) sin(wr), (12) 

Jo 

The first term describes the unitary part of the evolution 



The latter involves a summation over the various decay 
channels labeled by m with corresponding time-dependent 
decay rates A m (t) and arbitrary time-dependent system 

operators C m (t). In the simplest case, the rates A m as well J \ UJ ) w u 2 i w 2' v / 



being the noise and the dissipation kernels, respectively. 
In this paper we choose the Ohmic spectral density with 
a Lorcntz-Drude cutoff function, 



as the Hamiltonian H s and the operators C m are assumed 
to be time-independent, that is, it is the Markovian case. 
Note that, for arbitrary time-dependent operators H s (t) 
and C m (t), and for A m (t) > the generator of the mas- 
ter equation ([7]) is still in Lindblad form at each fixed 
time t, which may be considered as time-dependent quan- 
tum Markovian process. However, if one or several of the 
A rn (t) become temporarily negative, which expresses the 

presence of strong memory effects in the reduced system Entanglement measure quantifies how much entanglement 



where 70 is the frequency-independent damping constant 
and usually assumed to be 1. w is the frequency of the 
bath, and uj c is the high-frequency cutoff. The analytic 
expression for the coefficients 7(t)and A(t) are given in 
[33]. 

3.2 Measuring entanglement 



dynamics, the process is then said to be non-Markovian. is contained in a quantum state. The entanglement mea- 
For the system considered in Fig. 1, two Josephson junc- sure for a state is zero iff the state is separable, and the 
tion qubits coupled by the common current-biased Joseph- bigger is the entanglement measure, then more entangled 
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is the state. By axiomatic approach, (i) it is any nonnega- The general form of an "X" density matrix is as follows 
tive real function of a state which can not increase under 



^ a m + in 



b e + if 
e-if c 
\m — in d j 



(15) 



local operations and classical communication (LOCC) (so 

.... P 

called monotonicity) (57) : (ii) it is zero for separable states; 

(hi) and/or it satisfies normalization, asymptotic conti- 
nuity, and convexity. There are operational entanglement n , , , . , , , , n 
J ' J ^ ° such states arc general enough to include states such as 

measures such as distillable entanglement, distillable key T1T n i j j 

° ' J the Werner states, the maximally entangled mixed states 

and entanglement cost, as well as abstractly defined mea- /,,„,,n \ , , , , ., , • . 

° ' J (MhMSs) and the Bell states; and it also arises m a wide 

sures such as ones based on convex roof construction (e.g., r , • , ., « , , 

v ° ' variety of physical situations. A remarkable aspect of the 

concurrence and entanglement of formation) or based on a » v „ r ■ j ; , • .i . ,i , r . u 

° ' X form mixed states is that the time evolution of the 

distance from a set of separable states such as the relative , ,. / c \ i , j r. av. • i «v» *■ 

^ master equation (5) determined by tne initial X form 

entropy of entanglement R5] . One of the most famous mea- j j ^ i a- mi ■ i r 

1J o i— i 1S maintained during the evolution. Itus particular torm 

sures of entanglement is the Wootters' concurrence [55] t ,, , ., , . , , , . nl , , 

° 1 — ' of tne density matrix allows us to analytically express tne 

of two-qubit system. We will use it to study the entan- „__„ „ + „„ 

concui rcncc at time t as 



glement dynamics and obtain the entanglement transfers 
[59 56 55.. For a system described by a density matrix p, 



CZ(t)=2max{0,K 1 (t),K 2 (t)}, (16) 



the concurrence C(p) is where K x (t) = y/e 2 (t) + f 2 (t) - y/a(t)d(t), and K 2 {t) = 



v /m 2 (t)+n2(t)- v^(i). 

C(p) — max(0, ~ \ZA~2 — >/A~3 — \/Xi), (14) ^ ne non-Markovian master equation © is equivalent 

to a system of coupled differential equations, the first four 

where A x , A 2 , A 3 , and A 4 are the eigenvalues (with A x the of which descri be the time evolution of the populations, 

largest one) of the "spin-flipped" density operator £ — namc ly 

PK ® a v )P*K ® °f )> where P* denotes the complex h = - 2 {A(t) + 7 (i))o + (A{t) - 7 (i))(6 + c) 
conjugate of p and a y is the Pauli matrix. C ranges in 
magnitude from for a disentanglement state to 1 for a 
maximally entanglement state. 



b = (A(t) + 7 (t))a - 2A(t)b + (A(t) - j(t))d - 2Jf 
c = (A(t) + 7 (t))a - 22\(i)c + (A(t) - 7 (t))d + 2Jf 
d = (A(t) + 7 (t))(6 + c) - 2(4(i) - 7 (t))d, (17) 



The initial state chosen is the previous generated maxi- 
mal entangled state, 1^) = (\g}i\e)2 — i\e)i\g)2)/V2, which 
is an "X" form mixed state [60161162] that has non-zero 
elements only along the main diagonal and anti-diagonal. e = — 2A(t)e + (Eji — Ej^jf 



while the other equations describe the time evolution of 
the coherence: 
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Fig. 2. Color online) Non-Markovian entanglement concurrence C p (t) dynamics as a function of the thermal energy "fcsT" for 
the conditions \Eji — Ej%\ = and \Ej\ — Ej%\ = 4Ej, respectively. 



/ = -2A(t)f - (E, n - E J2 )e + J(b - c) 
m = -2A(t)m + {En + E J2 )n 
h = -2A(t)n - {E J1 + E. n )m. (18) 

In the following subsection we will bring to light the fea- 
tures characterizing the dynamics of superconducting en- 
tanglement. 

3.3 Numerical demonstration 

Using the same experimental parameters as |12j : Cj = 
4.3pF,/ = 13.3/i/l, l h = 0.9725/ , Cj/C ~ 0.1 and 
ujq = 2ir x 6GHz, which is chosen as the norm unit. The 
temperature is regarded as a key factor in a disentangle- 
ment process. Another reservoir parameter playing a key 
role in the entanglement dynamics is the ratio r = uj c /ujo 
between the reservoir cutoff frequency ui c and the system 
oscillator frequency By varying these two parameters 
ksT and r, the time evolution of the open system varies 
prominently for different cases. 



At first, let's consider two extreme cases, the strongest 
coupling (<P k = 0,^ ), and the decoupling (<P k = <P /2), 
which represent Ej k = ±2Ej and 0, respectively. In Fig. 
2, the time evolutions of the non-Markovian system con- 
currence for various values of temperature are plotted in 
the two extreme cases \Eji — Ej 2 \ — and \Eji — Ej 2 \ — 
AE°j. Here, we choose the ratio r — 0.1 and ksT ranges 
from to lOOwo. We can also compare the non-Markovian 
entanglement dynamics in the above two cases clearly. Fig. 
2(a) is the case of \Eji — Ej 2 \ = 0, which means both the 
two superconducting qubits have the same phase <£fc . The 
oscillation of the concurrence is also displayed in Fig. 2(a) 
besides the entanglement sudden death (ESD) and en- 
tanglement sudden birth (ESB). The concurrence decays 
with small amplitude in Fig. 2(b), \Ej\ - E J2 \ = AE°j . 
Both (a) and (b) show that the lower the temperature the 
more prominently the entanglement. The sudden death 
and birth behaviors shown in Fig. 2 are new features for 
physical dissipation and is induced by classical noise as 
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Fig. 3. (Color online) The evolutions of the concurrence C(t) of the superconducting qubits in a low temperature environment 
with ksT = 0.03c^o for (a) ratio r = w c /wo = 0.3, (b) r = 1, and (c) r — 10, respectively. In these figures, we mark the solid 
line as \Eji — Eji\ = 4JSj, and dashed line as \Ej\ — Ejq\ = 0. Figure (d): entanglement dynamics in a high temperature 
environment, ksT = ujq, magenta dashed line r — 0.3, blue dashed line r = 1, and red dashed line r = 10. 



well as quantum noise. The oscillating phenomenon em- 
bodies the non-Mar kovian memory effect in Eq. (6). The 
concurrence descended when A{t) — j(t) > and ascended 
when A(t) - j(t) < 0. 

In Fig. 3, we show entanglement dynamics for (a) r <^ 
1 ( magenta line), (b) r = 1 (blue line), and (c) r ^> 1 
(red line) in a low temperature reservoir (ksT — 0.03o;o), 
which implies that the ratio r = Ld c /^a between the reser- 
voir cutoff frequency w c and the system oscillator fre- 
quency ll>o plays a key role in the dynamics of the sys- 
tem. Obviously, the larger the ratio r the shorter the 
concurrence lasting time. Moreover, like Fig. 2, the con- 



currence decays with small amplitude in the condition 
\Eji — Ej2\ = 4Ej but has the same lasting time. The 
similar results can be found in [63164) . In Fig. 3(d), en- 
tanglement dynamics in high temperature environment, 
ksT — wo, magenta dotted line r = 0.3, blue dotted line 
r = 1, and red dotted line r = 10. In this figure we can 
see the ESD and ESB phenomenons obviously. From these 
simulations, we find that the entanglement can be open- 
loop controlled and the ESD time can be prolonged by 
adjusting temperature, r and the superconducting phases 
<Pk in the superconducting qubit systems. 



Wei Cui et al.: Non-Markovian entanglement dj 

4 Conclusions 

In the present work, we have theoretically studied the 
entanglement generation and dynamics in coupled super- 
conducting qubits systems. We characterize the entangle- 
ment by the thermal energy ksT, the ratio r and the en- 
ergy levels split of the superconducting circuit: \Eji(4>i) — 
^72(^2)!- Non-Markovian noise arising from the struc- 
tured environment or from strong coupling appears to 
be more fundamental and the ESD and ESB phenomena 
are analyzed in this paper. Our simulation results demon- 
strated that the lower the temperature the more promi- 
nent the entanglement. Moreover, the ESB phenomenon 
embodies the non-Mar kovian memory effect. We also find 
that the entanglement can be open-loop controlled and the 
ESD time can be prolonged by adjusting the temperature, 
r and the superconducting phases <Pk in the superconduct- 
ing qubit systems. Superconducting qubits offer evident 
advantages due to their scalability and controllability. We 
hope that such techniques will be experimentally imple- 
mented in the near future. 
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